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Abstract 

A number of exact results for two-loop three-point diagrams with massless in- 
ternal particles and arbitrary (off-shell) external momenta are presented. Divergent 
contributions are calculated in the framework of dimensional regularization. 
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1. There are several reasons why evaluation of loop diagrams with massless internal 
particles is important, namely: (i) some particles are really massless (photons, gluons); 
(ii) masses of some particles can be neglected in high-energy processes; (iii) such results 
are needed for constructing the asymptotic expansion of massive diagrams in the large- 
external-momenta limit; (iv) examination of such diagrams helps us to develop and apply 
new methods of loop calculations. 

There are two basic "topologies" of the three-point two-loop diagrams: planar (fig. la) 
and non-planar (fig. lb) ones. It should be noted that the diagram in fig. lb is symmetric 
with respect to all external lines, while the diagram in fig. la is symmetric with respect 
to the two lower lines only and therefore should be considered together with two other 
diagrams, corresponding to permutations of external momenta. In realistic calculations, 
however, we can also get some numerators (connected with spinor or vector structure 
of the particles). As a result, some of the denominators can cancel in separate terms 
corresponding to the decomposition in terms of the scalar integrals (some examples can 
be found in ref. jl|). Such integrals (that can be obtained from the diagrams of the types 
shown in fig. l(a,b) when some of denominators are removed) are shown in fig. l(c-f). All 
of them are shown for the case when they are symmetric with respect to the two lower lines, 
and the results for "rotated" diagrams can be obtained by permutations of the external 
momenta (as in the case of the diagram in fig. la). The diagrams in fig. l(d,e,f) are 
divergent, and we shall use dimensional regularization || (with the space-time dimension 
n = 4 — 2e). The main difficulty is that we need to evaluate arising one-loop triangle 
diagrams up to terms of order e. 
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Figure 1: 

The remainder of the paper is organized as follows. In section 2 we collect some 
useful formulae for one-loop triangles. In section 3 we consider planar diagrams (shown 
in fig. l(a,c)), while section 4 deals with the non-planar graph (fig. lb). Then, section 5 
examines the contributions of fig. l(d,e,f). In section 6 (conclusion) we discuss the results. 

2. In this section we shall briefly summarize some useful formulae for one- loop triangle 
diagrams. We shall employ them below, when evaluating all two-loop-order contributions. 
The scalar one-loop three-point Feynman integral with ingoing momenta Pi,P2,P3 
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(such that pi + P2 + P3 = 0) is defined, in n dimensions, as 

d n r 



J(n; z/i,z/ 2 ,z/ 3 | pl,pl,p 2 3 ) = J 



((p 2 - r) 2 )^^! + r ) 2 ) v z(r 2 y* 



The powers of propagators, z/j, can be non-integer if we employ analytic or dimensional 
regularization schemes (in the latter case we can obtain i/'s that depend on n as a result of 
previous loop integrations). Here and below the "causal" prescription {p 2 )~ u <-> (p 2 + i0)~ u 
is understood. Usually we shall omit the arguments p\-,p\-,p\ in J. 
The Feynman parametric representation of (|l|) is 



J(n; z/i,z/ 2 ,i/ 3 ) = 7r n/ i 



n/2;l-n r (E^-7l/ 2 ) 
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If the sum of the powers of denominators, v % = v \ + ^2 + v 3 , is equal to n, the triangle 
is "unique", and a very simple result can be obtained from (fj) (see ref. 0) 
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There are also two other useful relations of such type (see ref. H): 
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For 1/1 = z/ 2 = 1/3 = 1, n = 4 we define 
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where the dimensionless variables x and 7/ are 
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The function ^\x,y) can be represented as a Mellin-Barnes integral (see in [H^ 
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or as the parametric integral, 
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Figure 2: 



where the denominator can be also represented as a propagator, 

pi (y? + (l- x -y)£+ x ) = ( Pl + & 2 ) 2 
The result in terms of dilogarithms is: 

' py 



$W {x, y) = \ I 2 (Li 2 (-px) + Li 2 {-py)) + In V - In I 
A I x 1 + 



+ ln(px) ln(py) + 
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where 



-M^, y) = \J( l - % - v? - key , p(x,y) = 2(1 
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The results of the type of flTTD are well-known (see, e.g., in ref. ||). Here we use the 
notation of the papers |7], §, [|. It should be noted that the same function $W also occurs 
in the results for some other diagrams (for example, the one-loop box diagram M, |9[] and 
two-loop vacuum diagram with masses [0, 11]). Moreover, in four dimensions one-loop 



diagrams with larger number of external lines can be reduced to the four-point function 
(see, e.g., in refs. H12|), and therefore to a combination of $^ functions. 

3. In this section we will present the results for the diagrams shown in fig. l(a,c). The 
"ladder" diagram (fig. la), 



c^(pIpI,pI) 
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® {2 \x,y), 
(13) 

was calculated in ||. In this paper we will show that the result for the diagram in fig. lc 
can be expressed in terms of the same function. 

Let us introduce analytic regularization of these diagrams by shifting the powers of 
propagators by 6. t (as shown in fig. 2), provided that 61 + 62 + 63 = 0. Since these diagrams 
are convergent, the limit 6i — > should give results that do not depend on the way how 
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these <5's vanish. By use of the "uniqueness" relations ([3]) and @, the diagram on the 
r.h.s. of fig. 2 can be reduced to one- loop integrals ([[]) (see in [|J), and the result is 

in 2 „T(l-5i) f 1 



(p 



±-(pj) 5 ^ J(4; 1, 1, 1 - &) + -^-(p 2 )* 2 J(4; 1, 1, 1 - fc)} . (14) 
1O3 M3 J 
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Then, using the parametric representation 



^ 2 1 K t (vo-* -(*/*)- 



J(4; 1, 1, 1 + 5) = . ^t- x — / d£ -^f- (15) 







and considering the limit Si — ► (X) 5, = 0), we obtain |§ 

* (2,fe »> = ~\ I ^ + il -t-m + x lnf ( ln I + bf ) ( ln I + 2 lnf ) • (16) 

This integral can be evaluated as 

$ (2) (x, y) = j |6 (Li 4 (-px) + Li 4 {-py)) + 3 ln | (Li 3 (-px) - Li 3 {-py)) 

+\ In 2 ^ (Li 2 (-px) + Li 2 (-py)) + - ln 2 (px) ln 2 (py) 
2 x 4 

+ y ln(px) ln(py) + ^ In 2 | + , (17) 

where the polylogarithms Li^ (z) are defined as (see ref. jjTTJ) 

Ll Af (2) = 777 777 / d £ 
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Note that in ref. || Q these results fll6|) and (fPTD were generalized to the case of ladder 
diagrams with an arbitrary number of loops. Moreover, it was shown that the four-point 
ladder diagrams can be expressed in terms of the same functions as three-point ones (for 
example, the result for the "double box" diagram contains the same function $( 2 - ) ). 

Let us consider now the diagram on the l.h.s. of fig. 2 (it corresponds to the diagram 
in fig. lc). Using the relation @ for the lines that are connected with the "internal" 
vertex of this diagram, and applying also the formula (|3|), it is easy to obtain the result 
of the same form as (0), multiplied by the factor that is indicated in front of the r.h.s. 
of the equation in fig. 2. Considering the limit <5j — > (J2$i — 0) we find that the result 
for the diagram in fig. lc (at n = 4) is 

plC^(plplpl) = ^^(x,y), (19) 

P3 

where ^ is defined by (pD, (0). The equation ([H) (see fig. 2) can be easily checked 
by joining external lines with the momenta p\ and p 2 together and integrating over the 



5 



obtained loop momentum. The fact that both three-loop propagator-type diagrams should 
give the result proportional to 20((5) is well-known |l5j (the result for the ladder graph 
was calculated in [|1(], [IT]). It should be noted that the equation shown in fig. 2 can be 
generalized to the case of larger number of loops. 

4. In this section we shall examine the non-planar graph shown in fig. lb. The 
corresponding Feynman integral is defined as 



d n q d n r 



q 2 r 2 (p 1 - q) 2 (p 2 + r) 2 (p 1 - q - r) 2 (p 2 + q + r) 2 



(20) 



This function is symmetric with respect to all external lines, and each single loop of this 
diagram corresponds to a four-point function. 

To evaluate the crossed diagram , it is convenient to use Fourier transform to the 
coordinate space. It is easy to show that the "topology" of the diagram remains the same 
in the x-space (see fig. 3): there are three points X\, x 2 , x 3 (associated with external 




X2 



Figure 3: 



vertices), and there are two internal points z\ and z 2 , each of them being connected with 
all the three Xj. Note that in the x-space both integrations (with respect to the positions 
of Z\ and z 2 ) are independent. Therefore, the crossed diagram in fig. 3 factorizes (in the 
x-space), and the Fourier transform of (|20|) gives 



d 4 z 



(27T) 1 



Z 2 (X\ — X3 — z) 2 [x 2 — X3 — zf 



(21) 



Note that the integral in braces has the same structure as the momentum- space one-loop 
integral ([J). So, we may apply the Mellin-Barnes representation (H) to each of the two 
factorized integrals in (pT|). After we have done it, we get the following x-dependent 
structure in the integrand: 



(x 2 - x 3 ) 



u\+U2 r 



K X X - X 3 ) 



V1+V2 r 
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where u%, V\ and u 2 , v 2 are the Mellin-Barnes contour integral variables. In the momentum 
space, this product (|22|) corresponds to the "unique" triangle 



J(4; 2 — «i — u 2) 2 - Vx - v 2 , u% + u 2 + v l + v 2 ) 
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that can be written as a product of the powers of pf (see eq. (|3|)). Finally, collecting 
all gamma functions occurring in the transformations involved, we see that the resulting 
four-fold Mellin-Barnes integral (in the momentum space) can be split into two two-fold 
integrals, each of them corresponding to the one- loop function @! So, we obtain a very 
simple result: 



ITT 

,pl 



2\ 2 



$«(x,y) 



(23) 



with $W(x,y) defined by (|)-(|TT]) (therefore, can be expressed in terms of diloga- 
rithms and their products). We have mentioned this factorization property in ref. [18 . 
This result (^) has also been checked by joining two external lines and integrating over 
the loop, that gives the correct result proportional to 20^(5) (see ref. Hl6f). 



5. The diagrams shown in fig. l(d,e,f) are divergent, and we shall employ dimensional 
regularization |2[ to calculate the singular and finite (in e = (4 — n)/2) parts. All these 
diagrams contain one-loop propagator "bubbles" that can be easily calculated as 



i 1+2£ 7r 2 - £ T 2 (l-e) T(e) 



i 1+2e 7r 2 - £ 



(p 



2\e 



T(2 - 2e) 



(p 



2\e 



r(i + e) 



e 6 



(24) 



where p is the "external" momentum corresponding to the "bubble", and we put the 
dimensional regularization scale parameter p, = 1. The resulting integrals correspond 
to eq. (0), and they are: J(4 - 2e; 1, 1, 1) (fig. Id), J(4 - 2e; 1,1,1 + e) (fig. le) and 
J (4 — 2e; 1, 1, e) (fig. If). The problem is that these triangle integrals should be evaluated 
up to their e-parts (because they are multiplied by Q24D containing 1/e). 

We start with the diagram shown in fig. Id. We employ the "uniqueness" relation (^) 
to make the following transformations: 



J(4-2e; 1,1,1) 



r(i + e)r 3 (i-e) 



{pi)- e J(4-2e;l-e,l-e,l + e) 



r(l - 2e) 

= (plr £ (Pi)' 6 J(4-2e;l,l,l-2e). 



(25) 



Using these formulae and the Feynman parameters (fj), we can obtain the following one- 
parametric integral representation: 



J(4-2e; 1,1,1) 



7r 2 - £ i 1+2£ r(i + g )r 2 (i- £ ) i } der e ((2/Q- e -(s/Q- 

{plY+c T(l -2s) eJ (ye + {l-x-y)t + x 



t4- (26) 



From here, we immediately get a correct result ©, @) for the case of e = 0. Note that 
the denominator in fl26|) is the same as in (|) and (p!5|), and it can be represented as flTOD 
(this makes it possible to use (p6|) as a "block" in multiloop calculations). 
Let us denote 



J(4-2e; 1,1,1) 



7T 2 - £ i 1+2£ 



(pi) 



r(l + £ ) (x, y) + e (x, y) + 0{e 2 ) } 



(27) 
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where is defined by @, fllTI). The parametric representation for ^I^ 1 ** can be obtained 
from dffi), and it is 

* w (*. »> = 7 { ( in i + 2in «) Hye + (i - x ~ v)( + x) 

-21n?/ln£-21n 2 £- - ln(xy) In -} . (28) 

This integral can be evaluated in terms of polylogarithms (|18|) (up to the third order), 
and we get 

*<"(*. y) = -I ( 4 Li 3 (_^1±M) + 4 Li 3 (_MI±^) _ 4 Li 3 f f^) 
A [ \ 1 — /rxy y y 1 — p^ar?/ y y 1 — p z xj/ J 

+2 Li3 ( p^l) + 2U ( eaQ+M) _ 2 Li3 ( Ay ) _ 2<(3) 

V 1 + p x J \ 1 + py J v ' 

-2 In y Li 2 ( pX{1 + PV) ) - 2 In x Li 2 f Mi±^l ) 

2 2 2 

-- ln 3 (l - p 2 xy) + - ln 3 (l + px) + - ln 3 (l + py) + 21npln 2 (l - p 2 xy) 

-2\n(l- p 2 xy) ( hxipx) ln(py) +ln-ln 1+PV +21n(l+px) ln(l+py) + — J 
\ x 1+px 3 y 

+1 ln(p 2 xy) (ln(px) Hpy) + In V - In \±™ - In 2 + , (29) 

2 y xl + px 1 + px 3 y J 

where A and p are defined by flT2"|). This expression is explicitly symmetric in x and y. 
Note that in ref. |19| the massive case of such diagram was considered, and the result for 
the £-part also contained trilogarithms. By use of the integration-by-parts technique [J20| , 
we can obtain the results for the integrals J(4 — 2e; Ni, N 2 , N 3 ) (with integer N\, N 2 , N 3 ) 
in terms of the same functions @ (moreover, (|27|) is the only integral of such type where 
iffC 1 ' occurs in the e-part, other integrals will contain $ < - 1 - ) only). 

Another integral that we need to know up to the e-part, is J(4 — 2e\ 1, 1, 1 + e) (it 
corresponds to the diagram in fig. le). If we write Feynman parametric representation for 
this integral, we may observe, that the e-part of the parametric integral can be represented 
as a linear combination of e-parts of the integrals J(4 — 2e; 1, 1, 1) and J(4; 1, 1, 1 + e). 
The result for the latter integral follows from eq.(|15|) and was obtained in || (it can be 
expressed via the ^^-function). Finally, we get 

J(4-2e; 1, 1, 1+e) = ^^T(l+e) {d>«(z, y) (l - | ln(xy)) + y) + (9(e 2 )} , 

P3 (30) 
with the same as in (0). 

The result for J(4— 2e\ 1, 1, e) (corresponding to the diagram in fig. If) can be obtained 
either by direct evaluation (like (|26|)) or in a way similar to the evaluation of (0). It is 
divergent, but does not contain the function: 

7T 2 - £ i 1+2£ 1 

J( 4-2 ei i,i, e ) = - 5i55 -r(i + e ) 5?r -^ 
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|~ _e Inxlny + e (1-x-y) ® {1) {x,y) - e y + 0(£: 2 ) j 



(31) 



The results fl3"0|) and (|3T| ) enable one to evaluate other integrals J(4 — 2e; N\, N 2 , N 3 + e) 
(with integer N ly N 2 , N 3 ) by using the formulae of the paper |7j]. For example, 



7l 2 ~ £ i 1+2£ _ , 1 1 



J(4-2e;l J l,2 + e) = - r _ 5 -r(l + e 



(p|) 2+2e v 7 2(1 + e) xy 
x | h 21nx + 2\ny - e (21n 2 x + 2 In 2 ?/ + lnxlny) 



7T 2 1 

-3e (1 - x - y) (x,y)+s — + 0(s 2 )\ 



(32) 



It should be noted that the \ je pole in (|3l]) has an ultraviolet origin, while the 1/e pole 
in (|3^) has an infrared origin. 

6. In this paper, we have presented exact results for all the types of three-point off- 
shell contributions shown in fig. l(a-f). An interesting fact is that some diagrams can be 
expressed in terms of the functions corresponding to other diagrams. For example, the 
results for the diagrams shown in fig. l(a,c) involve the same function (see also fig. 2), 
and the result for the non-planar diagram in fig. lb can be represented as the one-loop 
triangle integral squared. There are only three independent functions occurring in the 
results (up to the two-loop level), namely: $ < - 1 - ) ([□]) (in fact, this function is connected 
with the one-loop triangle), $^ 2 ^ flTFD and All these functions are evaluated in 

terms of polylogarithms (]T8|) . 

The presented results are necessary for constructing asymptotic expansions of two-loop 
three-point massive diagrams in the limit of large external momenta (see, e.g., in [21] and 



references therein). The two-loop two-point case of such expansion has been considered 
in ref. ||11|| . They can also be used for evaluating some radiative corrections (like two- loop 
off-shell contribution to the three-gluon vertex). We should remember, however, that 
some scalar products (that may occur in the numerator) cannot be represented in terms 
of the denominators involved, and such terms should be considered separately. It should 
be noted that the on-shell singularities (that occur when some of the external momenta 
vanish) can be parametrized by putting pf = fi 2 (for those momenta that vanish), and 
they will appear as powers of ln/i (up to the fourth order). If dimensional regularization 
is used to regulate these on-shell singularities, they appear as poles in e up to 1/e 4 p2| . 
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